Disorder Effects in Superconducting Multiple Loop Quantum Interferometers 
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A theoretical study is presented on a number N of resistively shunted Josephson junctions connected 
in parallel as a disordered ID array by superconducting wiring in such a manner that there are N — 1 
individual SQUID loops with arbitrary shape formed. Under a constant current bias I > I c , and 
irrespective of the degree of the disorder, but depending on the strength of magnetic field B, all 
junctions in the array oscillate at the same frequency vb- Computer simulations of the full non 
linear dynamics of a disordered junction array reveal: (i) the frequency ub is not a periodic function 
of B, (ii) in the overdamped junction regime vb displays a sharp global minimum around B — 0. 
For zero inductive coupling the problem becomes equivalent to a virtual single junction model. 
PACS numbers: 85.25.Dq, 85.25.Am, 85.25.Cp 



So-called weak links, or Josephson junctions, are the 
basic active elements of superconductor quantum elec- 
tronics. A key feature of a weak link between two super- 
conductors, 1 and 2, is the property that there can flow 
a dissipationless macroscopic supercurrent I s (tp) due to 
the tunneling of Cooper pairsu with charge 2e. This su- 
percurrent depends on the gauge invariant phase differ- 
ence tp = 61 — 62 + |f J-, 2 < ds, A > of the macroscopic 
BCS pairing wavefunctions on either side of the weak 
link. Joseohson junctions made with modern fabrication 
techniquesa often have a sandwich type layered geome- 
try, with a thin non superconducting tunneling barrier 
in the middle between two thick superconducting elec- 
trodes. In recent time also other types of weak links, 
for example of the bi-crystal type,rJpecame important in 
high-temperature superconductorsaa. For an ideal S-I-S 
junction the supercurrent is connected to the phase differ- 
ence tp across the tunneling barrier by I s (tp) — I c sin (p. It 
is important to realize that the supercurrent 7 S flows sta- 
tionary provided it does not exceed a characteristic crit- 
ical current I c , the so-called Josephson critical current, 
which determines the maximum dissipationless current 
that can flow across a tunneling barrier. In general, I c 
depends on the material properties of the junction, on 
temperature T, and also on magnetic field B = rot A. 
Applying to a Josephson junction a bias current I with 
a constant strength I > I c , there appears a rapidly oscil- 
lating voltage signal V (t) across the junction, which de- 
termines the rate of change of the time dependent phase 
difference tp (t) according to 



hd t tp(t) = 2eV(t) 



(1) 



This is the fundamental non stationary Joseph- 
son relation which governs the physics of weak 
superconductivity!! So, for I > I c there flows, besides 
the dissipationless supercurrent I s , also a dissipative nor- 
mal current in the junction, whose physical origin is 
the transfer of single (unpaired) electrons. _ 

Within the range of validity of the RCSJ modelO, the 
dissipative current may be described with sufficient ac- 
curacy as a superposition of an ohmic current, charac- 
terized by a parallel ohmic shunt resistance R, and a 
displacement current, which is characterized by a paral- 



lel geometric shunting capacitance C describing electric 
polarization inside the tunnelling barrier. The total junc- 
tion current I is then: 



I = I c sin 99 (t) 
The time average 
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is the dc voltage part of the in general not sinusoidal 
voltage signal V(t) across the electrodes of a Josephson 
junction. For a strongly overdamped junction, C = 0, 
one finds, assuming a constant bias current I > I c , a 
relatively simple formulaQ: (V) — R JP — . The dc 
voltage (V) is connected to the oscillation frequency uj — 
2-kv of the voltage signal V (t) by: 



hp = 2e (V) 



(4) 



This result for the voltage response function (V) of a 
weak link suggests a spectroscopic interpretation. When 
a Cooper pair is transferred from the superconducting 
side 1 to the superconducting side 2 of the junction, un- 
der conditions where I > I c , a microwave photon with 
energy 2e (V) is released in the form of one quantum of 
electromagnetic radiation (Josephson radiationtl) . 

As far as macroscopic quantum interference is con- 
cerned, it was actually knowru long before the discovery 
of the Josephson effects, that magnetic flux threading 
the area of a superconducting ring, made out of a mate- 
rial that is thick compared to the magnetic penetration 
depth A, should be quantized in units of the flux quantum 
$n = — 

2e ; 

Technical applications of the physics of weak 
superconductivityHB include ultrasensitive quantum in- 
terferometers, which indeed combine the afore mentioned 
Josephson effects with flux quantization. Consider, as 
indicated schematically in Fig. (la), a standard two junc- 
tion SQUID (for simplicity with symmetric junction pa- 
rameters) under the dc current bias I > 2/ c ._ Such 
a device is actually a flux-to-voltage transducer^. Let 
$ = (B, a) = |B| |a| cos a be the magnetic flux thread- 
ing the orientated area element a of the superconducting 
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SQUID loop, where a is the angle between the normal 
vector of the orientated area element and the magnetic 
field vector B, as depicted schematically in Fig. (la). The 
total magnetic field, B = B^+B 1 - 2 -', is then a superpo- 
sition of the primary external magnetic field B^ 1 ) , which 
generates the flux $W = (B^, a) one wants to detect, 
and a secondary magnetic field B^ 2 ) that results, for 
example, from the inductance L (or other impedance ef- 
fects) in the circuit. The screening current I sc circulating 
in the SQUID loop leads to a total flux $ = $W + $ (2 ). 
Dependent on the secondary flux term, $' 2 ) = — L -r-Jsc, 
there exists an optimal size \sll \ for any SQUID loopB. A 
dimensionless measure for the inductance of such a loop is 
0L = Note also that a two junction SQUID cannot 
be directly employed as a detector of absolute strength of 
external magnetic field. This is because the voltage re- 
sponse function (V xy ) of the SQUID, i.e. the time average 
of the rapidly oscillating voltage signal V xy (t) across the 
nodes x and y of the circuit, is a periodic function of the 
strength of external magnetic field, see Fig. (la). 

A straightforward extension of the standard two junc- 
tion SQUID is sketched in Fig. (lb.). This is a ID array 
of TV adjacent Josephson junctional connected in parallel. 
In particular, the area elements of the N— 1 SQ UID loops 
formed in this manner are all equal, e.g. a n = ell for all 
n. The voltage response signal (V xy ) vs. strength IB^I 
of external magnetic field of this array has the same pe- 
riod than a standard two junction SQUID with loop area 
|a L |, see Fig.(lb). 

A more general quantum interference device is ob- 
tained when the area elements a„ of the N — 1 loops 
in the array differ in size and, possibly, in orientation, 
as depicted schematically in Fig.(lc). If the sizes |a„| 
of the orientated area elements a„ of the individual su- 
perconducting loops are chosen in a disordered fashion 
the voltage response function (V xy ) vs. I B^ 1 ) I becomes 
nonperiodic, see Fig.(lc). Taking into account inductive 
couplings among the currents in the circuit, the maxi- 
mum loop size in the disordered array should coincide 
with the corresponding optimal loop size of a standard 
two junction SQUID, i.e. max|a„| = |ax,|. The voltage 
response signal (V xy ) vs. strength of magnetic field of a 
disordered junction array is, under a suitable dc current 
bias J, indeed a unique function of |B| around its narrow 
global minimum at |B| = 0. This suggests that it should 
be possible, e.g. by measuring control current (s) flowing 
through the wires of a set of suitably orientated compen- 
sation coil(s), to reconstruct absolute strength, orienta- 
tion and even the phase of an incident primary magnetic 
field signal, i.e. to determine the full vector BW(t). 

The n-th Josephson junction in the array has, within 
the range of validity of the RCSJ model, optional indi- 
vidual junction parameters C„ and I c ,n- The cor- 
responding current /„ flowing through the n-th Joseph- 
son junction is, according to Eq.(|^), determined by the 
gauge invariant phase difference (p n (t) across that junc- 
tion. The total current / flowing through the nodes x 
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FIG. 1. Voltage response (V xy ) in units of I c Rvs. external 
flux through largest area element a_L of interferometer for 
bias current I — 1.1 NI C and inductance fit = 0: a) symmet- 
rical SQUID (N = 2), b) periodic ID array (N = 11), c) 
disordered ID array (N=18), but with same total area as b). 

and y, respectively, of the circuit is then obtained from 
Kirchhoff's rule as the phase sensitive superposition of 
the individual junction currents /„: 



N-l r 



n=0 



I C) „ sin </>„(*) + (— — d t + TTS-Ot ) Vn (t) 



(5) 



Note that the gauge invariant phase differences ip n of ad- 
jacent Josephson junctions in the array are not indepen- 
dent, but are connected to each other by the condition 
of flux quantization: 



27T 
$0 



(B , a„) mod 27r 



(6) 



Here |a„| is the area of the superconducting loop con- 
necting adjacent Josephson junctions numbered as n and 
n — 1, respectively, and B denotes the magnetic field 
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threading the orientated area element a„ of this loop. 
Note that Eq.(||) applies quite generally, provided the 
superconducting material, out of which the connecting 
loops are made, is thick compared to the magnetic pen- 
etration depth A. In this case there exists a path inside 
the wire connecting, say, junction n with its neighbor 
junction n — 1, on which the superfluid velocity field v s 
becomes negligible small. So, TiWQ — along this 

path. Since all junctions in the array are connected in 
parallel, the rapidly oscillating voltage V n (t) at the elec- 
trodes of a particular Josephson junction, numbered as 
n in the array, is related to the signal V xy (t) between the 
nodes x and y of the circuit by 



responsive function of strength and orientation of mag- 
netic field, and it is strongly affected by the choice of the 
individual area elements a m . In general \Sn (B)| is very 
sensitive to permutations among the a m 's. 

In the overdamped junction regime, C = 0, under con- 
ditions where a constant current / is biased such that 
1 > \Sn (B)| / J/v = sin as, and assuming for simplicity 
a homogeneous static magnetic field B (and also time in- 
dependent area elements a m ), one finds an exact solution 
for the phase difference (po(t)'- 



V (t) = -d m (t) 



I C R 



j2 

N 



\S N (B) 



Jn + \Sn (B)| sin (lu b ■ t - a B ) 



V xy (t) = V n (t) + 



(da, E(i)) 



(7) 



x — >n — >y 



By Faraday's law the electric field E along an integra- 
tion path x — > n — > y, that starts at node x, traverses 
the tunneling barrier of the n-th Josephson junction just 
once, and then terminates at node y, is directly connected 
to the time derivative of the flux threading the area ele- 
ments of the ID array. Once the signal Vo(t) — j-dt(fo(t) 
is known the other voltage signals V n (t) across the elec- 
trodes of the n-th junction follow from 
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V n (t)-V n -i(t) = -d t (B(t) ,a„) 



(8) 



Takingpinto account the Biot-Savart type inductive 
couplingsB among the currents flowing in the circuit pro- 
hibits further simplification. However, it follows directly 
from Eq. (||) that one can indeed eliminate from Eq. (||) all 
phase variables ip n (t) in favor of a single phase, for ex- 
ample tpo(t), provided the extremely simplifying assump- 
tion is made, that all inductive couplings vanish. In this 
case the problem of N coupled Josephson junctions is 
mapped onto a virtual single Josephson junction model. 
With i = lEto 1 jfc, T N = A_2_ an d (f) = fo( i) 
there results a scalar differential equation determining 
the phase difference </> (t) : 

\S N (B)\ sm^(t)+6 N (B)}+T N (RCd 2 + d t )cf > (t) ( 9 ) 
= J N -^T N (RCd? (B(i) ,a c )+ d t (B (t) ,a R )) 



where we have defined (ao = 0) 
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N-l 



1 X -s R X ^ 

N RZ ^ 



2ni 



E 

m=0 



B , a rl 



(10) 



m=0 



n=0 



C ^ 

m— 



and J N = jjj-, I c = ]y 2n=0 I c,n, C = jj J2n=0 Cn- 

The complex function S N (B) = \S N (B)| e l5jv(B) denotes 
the characteristic structure factor of the ID Josephson 
junction array, as defined in Eq.(|lO|). It is an extremely 



For a static magnetic field B the voltage response func- 
tion (V X y) measured between the nodes x and y of the 
circuit is equal to the dc-part of the rapidly oscillating 
voltage signal Vo(t). All Josephson junctions in the ID 
array oscillate at the same frequency = 2-kvb, which 
is related to (V xy ) by: 



2e 



v B = (V ) 



IcR\U 2 N 



(B)| = (Vxy 



(11) 



Note that the oscillation frequency vb of such a local 
oscillator is even more sensitive to changes of strength 
or orientation of the external magnetic field B than the 
structure factor of the array itself, since |5jv (B)| enters 
Eq.([lll) quadratically. 

Consider, as a special case, an ordered array, consist- 
ing of N— 1 identical SQUID loops, such that (B , a„) = 
$ = (B , a^), and I c n = I c independent on the junction 

index n. Then the structure factor Sn (B) = sffl be- 
comes a simple geometrical series: 



N — 



N sin ( 7r 



exp 



m-(N-l) 
<P 



(12) 



In Fig. (lb) one observes the usual narrowing propor- 
tional to of the width of the voltage response signal 
(Vxy) around its minima. Note the periodicity property 
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N 



for all N > 2. For N = 2 Eq.(|ll|) is 

the periodic voltage response of a symmetric-two junction 
SQUID in the overdamped junction regimes. 

A structure factor with a much longer period may be 
obtained in a parallel junction array where the orientated 
area elements increase in size according to a linear rela- 
tion: 



a m = (2m - 1) ai 



(13) 



For simplicity, identical junction parameters R n , C n and 

Ir 



L Ci „ are assumed. Then: 
N-i 



Sn ^ = jf E ex p 



n=0 



(14) 
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The total area occupied by such a Gaussian array is 
(N — l) 2 ai, where ai is the smallest area element, and 
aAT_x = (2N — 3)ai is the largest area element. Com- 
pare a Gaussian array, with N — 1 area elements as de- 
scribed in Eq. (|l3|) , with a periodic array, consisting of 
Np — 1 identical SQ UID- loops with size For a use- 
ful comparison, both arrays should occupy the same to- 
tal area: (N — l) 2 ai = (Np — l)a^. Also the largest 
area element in the Gaussian array should coincide with 
the area element of an optimal single SQUID-loop, i.e. 
ajv-i = sll- Both requirements together imply for 
Np 3> 1 that the Gaussian array has the double number 
of junctions compared to a corresponding periodic junc- 
tion array: N ~ 2Np. To determine the period of the 
Gaussian array consider a case where the flux threading 
the area of the smallest element, ai, is equal to a rational 
multiple of half a flux quantum: (B , ai) = 77-%- ■ Then 
the largest area element in the array, ajv-i, is threaded 
by a flux $m = (l — ^§7) M in this case the struc- 



ture factor Sjv(B) ^- 
result of C.F. GaussE 



may be determined using a 
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Note the periodicity |s<**'+* a "> 



■»(**r) 



J , with period 

$2N = (2iV - 3) $ - Remarkably, for M = 2, and N = 
N1N2 being the product of two prime numbers Ni and 
N2, there holds the factorization: 
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(16) 



Apparently, such Gaussian junction arrays are governed 
by the laws of number theory (quadratic residuestLH) . 

The long periodicity of the structure factor vs. flux 
$W threading the largest area element &l of Gaussian 
junction arrays is also visible in the calculated voltage 
response function (V xy ), irrespective of the degree of 
the inductive coupling represented by the parameter fti,. 
This is illustrated in Fig. (2). Note the asymmetry of 
(V xy ) under <I> — * — $ (for a constant bias current I) 
for finite inductive coupling. As far as disorder is con- 
cerned, we also find that (V xy ) in Gaussian junction ar- 
rays is very responsive to adding small random fluctu- 
ations to the size distribution of the area elements, so 
that (V xy ) becomes non periodic with a pronounced an- 
tipeak only around $ = 0. As /3l increases, the difference 
max (V X y) — min (V xy ) decreases, and the linewidth of the 
global minimum (V xy ) around $ = ceases in this case to 
scale proportional to . Note all this applies in the over- 
damped junction regime. Our results for weak damping 
will be published elsewherdlj. 

We hope that an experimental verification of the pre- 
dicted magnetic field dependence of the voltage response 
function of disordered ID parallel Josephson junction 
arrays will stimulate the development of new types of 
robust superconducting quantum interferometers, which 



would allow (for the first time) a technically rather simple 
precision measurement of absolute strength of external 
magnetic fields. 
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FIG. 2. Voltage response {V xy ) in units of I C R vs. exter- 
nal flux <E >< ' 1 - ) through largest area element a_L for a Gaussian 
array with TV = 18 (overdamped) junctions for bias current 
/ = 1.1 N I c and various inductive couplings: a) /3l = 0, b) 
Pl = 0.3 and c) /3 L = 0.7. 
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